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Abstract 

The formulas directly connecting parton distribution functions (PDFs) and fragmen- 
tation functions (FFs) at the next to leading order (NLO) QCD with the same quantities 
at the leading order (LO) are derived. These formulas are universal, i.e. have the same 
fN| ■ form for all kinds of PDFs and FFs, differing only in the respective splitting functions 

$h ■ entering there. 
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The extraction of PDFs and FFs from the experimental data is one of the important tasks 
of the modern hadron physics. The most simple and transparent way to do it is the QCD 
analysis of the data on measured asymmetries and cross-sections in LO QCD. One of the main 
advantages of such analysis is that the central values and uncertainties of measured asymmetries 
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D ■ and cross- sect ions directly propagate to the central values and errors of PDFs and FFs extracted 
r^i l from these data in LO QCD (see, for instance, Fig. 3 in pLj). At the same time the situation 
\ with NLO analysis is much more difficult because instead of simple algebraic equations (see, 
^ ; for example, Eq. (2) in Ref. PQ) one deals there with complex integral equations (like, for 
■ instance, Eqs. (10)-(13) in Ref. [2]) for finding PDFs (FFs) we are interested in. The standard 
q ■ way to solve this problem is to apply the QCD analysis based on the fitting procedure (see 
\Q '. [2] and references therein). However, there are unavoidable ambiguities inherent in a fitting 
! procedure which become especially important when the quality of the fitted data is rather 
^ [ bad (small number of points with large errors). These are arbitrariness in the choice of the 
_ ■ functional form (with a lot of varied parameters) of the fitted PDFs and FFs at initial scale 
and also ambiguities in the error band calculation (ambiguities related to the deviation of x 2 
profile from the quadratic parabola and to the choice of A% 2 determining the uncertainty size 
- see discussion on this subject in Ref. [2]). Thus, it seems to be very useful if one could obtain 
NLO (NNLO, ...) results on PDFs/FFs using the respective LO results as an input, without 
loosing, thereby, all advantages of LO analysis. 

We start with some necessary notation and definitions. For the flavor non-singlet and singlet 
quantities we introduce the notation Q NS and V = (Q s , G), where Q^s can be either q NS (non- 
singlet combinations of quark densities), or Aq^s (non-singlet combinations of helicity PDFs), 
or combinations of transversity PDFs Axqiq) = hi q> g,. . ., or D^ s ("non-singlet" combination 
of FFs Dg),. . ., while Qs can be qs, Aqs, Dg,. . ., G can be g, Ag, D^,. . .. In this notation the 
DGLAP evolution equations (see j3] for review) look as 

Q 2 dV(Q 2 ,x)/alQ 2 = (a s /2n)[p(°\x) + (a s /27t)P^(x) + 0(a 2 s )} ® V(Q 2 , x), (1) 
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where the convolution (<g>) is given by 



(A g) 



f 1 f 1 f 1 dy x 

B)(x) = / dx\ I dx 2 S(x — XiX2)A(xi)B(x 2 ) = / — A(—)B(y), 
Jo Jo Jx y V 



and analogously for Qjvs with the replacement P(x,a s ) —> P(x,a s ) = P <yQ \x) + (a s /2'K)P^(x) + 
0(a 2 ). Here P is 2 x 2 matrix with the elements P qq , P qg , P gg , P gg , and the splitting functions 
for unpolarized PDFs and helicity PDFs can be found in the review [I], for transversity PDFs 
- in the review [5] , for FFs - in Ref . [6] and references therein. 

Following [7] it is convenient to define the evolution operators E and E (2 x 2 matrix with 



the elements E qq , E qg , E gq , E gg ) as 



Qns(Q 2 , x) = E(Q 2 ,x) ® Q NS (Q 2 , x), V(Q 2 , x) = E(Q 2 , x) ® V(Q 2 , x) 
Here we are interested in the initial conditional 

E(Q 2 = Ql, x) = 5(1 - x), E(Q 2 = Q 2 , x) = 1 5(1 - x), 



(2) 



(3) 



which allow to evolve Qns and V from the initial scale Qq to an arbitrary scale Q 2 . 

It is also convenient to use, following [7], the evolution variable t = (2//3 ) In (a s (Ql) / a s (Q 2 )) 
instead of the standard variable ln(Q 2 //i 2 ). Besides, we introduce the notation 



A \lo= A ' 



A 



NLO 



= A 



(4) 



for any quantity A at LO and NLO, respectively. 

From now on we consider only the nontrivial singlet case. Transition to the simple non- 
singlet case will be easily done in the end of calculations by making the replacement of the 
matrices with the respective commuting quantities. 

In terms of quantities t and E the DGLAP equations are rewritten in LO as 



-n&(i,x) = P (0) ®E(t,x), 



while in NLO they look as 



d 
dt 



where 



^E(t,x)= P (0) (x) + — R(x) + 0(a 2 s 



R(x) = P (1) (x) - 4- p(0) M 



E(t,x), 



(5) 

(6) 
(7) 



Solution of (j3J) with the initial condition (j3J) E(t = 0, x) — 1 5(1 — x) reads [7] 



E(t,x) = Exp(P (0) (a;)£)= 1 5(1 - x) + t P (0) (x) + - P (0) (x) ® P (0) (a;) + . . . , (8) 

while to solve NLO equation flS]) one can apply the elegant method of Ref. [7] based on 
the analogy with the perturbative quantum mechanics (see Eqs. (5.47)-(5.54) in Ref. [7]). 



2 We do not consider the asymptotic conditions [7] E (E) — > E (E) as Q 2 — > oo (see Eq. (5.57) in 0), since 
we deal only with particular realization @ of the general conditions given by Eqs. (5.18) in [7]. 



2 



Operating in this way one obtains the general solution of ([6]) in the form (for a moment we 
omit x dependence and 5(1 — x)) 




Jt' 



Putting t — > oo in (j^J) one reproduces the solution (Eq. (5.54) in Ref. [7J) satisfying the 
boundary condition E — > E as t — > oo. In turn, putting t = in ([9]) one gets the solution 




satisfying the boundary condition (|3]) we deal with. 

The key point to proceed is the condition that all PDFs and FFs should take the same 
values in LO and NLO (as well as in NNLO,. . .) as Q 2 — > oo: 



Though this asymptotic condition seems to be intuitively clear, let us argue it in some detail 
because of its great importance for what follows. 

Imagine that two researchers analyse in LO (the first) and NLO (the second) the same 
"ideal" data - the data available with tremendous statistics even in the Bjorken "sub-limit" (so 
high Q 2 values are accessible that the Bjorken scaling violation becomes invisible even within 
extremely small uncertainties on measured asymmetries and cross-sections). For determinacy 
and simplicity let us suppose that they analyse the imaginary "ideal" polarized SIDIS data 
on pion production and extract the valence helicity PDFs Auy, Ady from the proton and 
deuteron difference asymmetries (see Ref. [8] and references therein) measured in the Bjorken 
"sub-limit". The first uses LO formulas Af~ n ~ ~ (4Au v - Ad v )/{4u v - d v ) and A^~^ ~ 
(Aiiy + Adv)/{uv + dy) (i.e., performs the analysis analogous to one of COMPASS [9]), and 
the second their NLO generalization (Eqs. (6-10) in Ref. [5]). Besides, for self-consistency, 
both imaginary researches do not use the existing parametrizations on uy,dy but extract these 
quantities themselves (as well as the integrated over cut in z differences D\ — D 2 of favored 
and unfavored pion FFs) using the same SIDIS data on pion production averaged over spin and 
studying the quantities F^ d3He s — F^ d3ne y where in both LO and NLO only uy, dy and 
D\ — D 2 survive. It is obvious that all terms with convolutions £§> (see Eqs. (6-10) in Ref. [8J) 
distinguishing NLO and LO equations for finding Auy, Ady and uy, dy, Di — D 2 just disappear 
as one approaches the Bjorken limit, so that comparing the results on these quantities obtained 
in the Bjorken "sub-limit" both researchers could not discriminate between them. 

So, let us pass to limit Ql — >■ 00 in Eq. fl2]) using the asymptotic condition ffTT|) . Then, on 
the one hand (NLO evolution) 

V(Q 2 , x) = E(t -00, x) ® V(Ql -t 00, x) = E(t -»■ -00, x) ® V(Qo ->> 00, x), (12) 

and, on the other hand (inverse LO evolution) 



3 On simultaneous determination of valence PDFs and Di — Di from the SIDIS data see, for example, jlO) . 



Qns{Q 2 -> 00, x) = Q, NS {Q 2 y 00, x) , V(Q 2 -> 00, x) = V(Q 2 -> 00, x). 



(11) 



V(Ql -> 00, x) = E(f -> 00, x) <g> V(Q 2 , x). 



(13) 
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Combining Eqs. ( TT2|) and ( TT3|) one obtains 



V(Q 2 ,x) 



lim E(t,z) <g> E(-i,x) 



V(Q 2 ,x). 



(14) 



Using Eqs. ([H]), (HOI) and the relation \im.Q2^ 00 (a s / a s ) = 1 we arrive at the connection formula 
between NLO and LO flavour singlet PDFs (FFs) V and V at the same finite Q 2 value 



V(Q 2 , x) = \l 5(1 -x)- fl^ dr e A)T / 2 E(r, x) ® R{x) ® E(-r, x) 

® Exp (-f Q In fg{ P(°)(x)) is V(Q 2 , x), 



(15) 



where all dependence on the unreachable infinite point Ql just cancels out. 

In the non-singlet case the relation ()15p is significantly simplified. The terms E(t, x) 
Exp(rP(°)(x)) and E( — r, x) cancel out each other in the integrand and one easily obtains 



Q NS (Q 2 ,x) = [6(1 -x) + ^ (|P (( %) - lP^( 
® Exp (- f Q In P(°) (x)) ® Qa^Q 2 , s). 



(16) 



Eqs. (|15|) and (fT5"|) connecting flavour singlet and non-singlet quantities in NLO with the 
same quantities in LO is the main result of the paper. Let us briefly discuss their practical use. 

There are not any problems with application of Eq. (|T6|) and the task of reconstruction of 
NLO non-singlet quantities from LO ones is reduced just to the trivial calculation of the integrals 
entering the convolutions ®. Indeed, the parameter e = — (2//3o) ln(a s /& s ) is very small even at 
the minimal (the lower boundary of the experimental cut on Q 2 is usually about 1 GeV 2 ) really 
available Q 2 values, so that one can achieve very good accuracy keeping only few first terms in 
the expansion Exp (eP^(x)) = 5(1 - x) + e P^\x) + (e 2 /2!) P^(x) <g> + . . . Certainly, 

the same statement holds for term Exp (eP'°'(i)) in Eq. f[T5"j) . but there arises an additional 
problem how to deal with the integral over r. As usual, the problem is easily solved in the space 
of Mellin moments. Notice that Q 2 independent integral over r in Eq. f|T5il just coincide^] with 
the quantity —U(x) in Ref. [7] (see Eq. (5.45) in [7j), which enters the solution of DGLAP with 
the boundary conditions limQ2_ ) . 00 E (E) = E (E) (see footnote 2). Then, applying the inverse 
Mellin transformation, one easily obtains instead of fll5p the formula suitable^ for numerical 
calculations 



V(Q 2 ,x) 



15(1 -x) + 



a s (Q 2 
2tt 



C+ioo 



X 



C—ico 



dn U (n) 

2iri v ; 



Exp(e(g 2 )p(°)(x))®V(g 2 ,x), (17) 



where 2x2 matrix U(n) = f^ dxx n ~ l U(x) is given by Eq. (5.41) in Ref. [7]. 

In summary, the formulas allowing to transform LO parton distribution and fragmentation 
functions to NLO ones are derived. To obtain these formulas we use as an input only the 
DGLAP evolution equations and the asymptotic condition that PDFs (FFs) at different QCD 
orders become the same in the Bjorken limit. Due to universality of this input the connection 



4 Using Eq. (5.28) in Ref. [7] for U and the obvious relation Q 2 d[ Exp ((2//3 ) lnd s P< >) ® V]/dQ 2 = 
one can immediately check that r.h.s. of Eq. (fT5"]) indeed satisfies the NLO DGLAP equation (JTJ) . 

In Ref. one can find the efficient algorithm for the numerical calculation of the integral over n (proper 
choice of the integration contour, etc. - see discussion around Eq. (3.2) in Ref. 
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formulas are also universal, i.e. they are valid for any kind of PDFs (FFs) we deal with. Besides, 
it is obvious that operating in the same way one can also establish the connection of PDFs (FFs) 
at LO (as well as at NLO) with these quantities at any higher QCD order (NNLO, NNNLO, 
. . .), and the only restriction here is the knowledge of the respective splitting functions. 
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